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Abstract 

Inspired by the latest results of ATLAS and CMS on the search for the standard 
model (SM) Higgs scalar, we discuss in this article the correlations between Higgs-boson 
properties, low-energy observables, such as B -^s7) Bg — ji^ iJi~ , and (5 — 2)^, and the 
dark matter (DM) relic density. We focus on the corners of the MSSM parameter space 
where the pp — t- /i — t- 77 signal is enhanced due to the presence of a light stau state. 
In this region tan/?. Ma, At, and ji take large values, and we find striking correlations 
between many of the considered observables. In particular, the B Xg^ branching 
fraction is enhanced, while the Bg — )• ji'^ rate tends to be below the SM expectation. 
In contrast, the Higgs-boson couplings show good overall agreement with the preliminary 
experimental determinations, the DM abundance is consistent with observation, and the 
discrepancy in {g — 2)^ is reduced. The predicted deviations and found correlations 
could be tested in the near future and hence may become very valuable as guidelines 
and consistency checks. 



1 Introduction 



Recently both the ATLAS |1| and the CMS [2] collaborations have announced the existence 
of a new bosonic state with a mass of around 125 GeV. This discovery defines a turning 
point in the history of elementary-particle physics: the almost 50 year-long hunt for the Higgs 
boson has come to a dazzling end. One of the most central problems in high-energy physics is 
now whether the properties of the newly observed resonance agree with that of the standard 
model (SM) Higgs scalar. In fact, a major experimental effort is directed towards shedding 
light on this question by measuring the various decay rates of the new particle as accurately 
as possible. While the achieved precision is so far insufficient to draw any final conclusion, one 
cannot help but noticing that the LHC data match well the SM prediction with one possible 
exception: the central value of the measured relative strength of the diphoton signal is too 
large by about 70% [l]|2]. If this enhancement survives further experimental scrutiny it may 
become the first convincing evidence of physics beyond the SM. 

An economic possibility to modify the effective coupling between the Higgs boson and two 
photons without altering its main decay modes is provided by extra vector-like leptonic states 
(see for example [3}jl5]). In the context of the minimal supersymmetric SM (MSSM) a light, 
maximally mixed stau can be shown to increase the loop-induced /;,77 coupling [s], hence 
leading to the desired effect. 



In this article we dissect the "light stau scenario" extending previous analyses 3l4 6 ,10 ,12 



16 that are similar in spirit. In particular, we will show that enhancements of the — )• — )■ 77 



signal associated to virtual stau exchange occur only in a narrow sliver in MSSM parameter 
space that features large (and positive) values of tan/3, the trilinear couphng At, and the 
H parameter. In this corner of phase space we investigate in detail the indirect constraints 
arising from the most relevant 5-physics observables, the muon anomalous magnetic moment 
{{g — 2)fj), and dark matter (DM). We observe striking and testable correlations between the 
individual observables. Our main findings are as follows: i) An enhancement of the diphoton 
signal of the wanted order necessarily leads to an increase of the B — )■ Xsj branching ratio 
of around 30%. With improved theoretical and experimental determinations of the inclusive 
radiative B decay this would represent an unambiguous signature of the light stau scenario. 
ii) For positive gluino mass parameters M3, the decay rate of Bg — )■ /i'''/^" tends to be below 
its SM prediction with the exact value depending sensitively on At and Future precision 
measurements of the purely leptonic Bg decay by ATLAS, CMS, and LHCb hence represent 
sensitive probes of the MSSM parameters and its particle spectrum. Hi) The supersymmetric 
(SUSY) parameters that play an important role in in the i?-physics observables also leave 
an imprint in the h ^ bb decay rate which typically turns out to be larger than expected 
from a SM Higgs boson. Tentatively, this seems to be in accordance with the findings re- 



ported by CDF 17 and D0 ^|. The h — t- r+r^ channel, on the other hand, is in general 
suppressed, which is again in line with observation [2]. iv) Under the assumption of a light 
slepton/sneutrino spectrum the predicted MSSM corrections in {g — 2)^ allow for a good de- 
scription of the experimental data, v) In the presence of a light bino-like neutralino the specific 
MSSM scenario under study is compatible with thermal DM. In fact, the interplay between 
an enhancement in /i — )■ 77 and the requirements imposed by the relic density essentially fixes 
the mass splitting between the lightest neutralino and stau. 
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The intriguing correlations found in our article can be used as guidelines for the direct 
searches of SUSY particles, and will become extremely valuable as consistency and cross 
checks, in case the LHC high-pr experiments will start to see the first scalar partners. 

This work is divided into two parts. In Section [2] we pursue an analytic approach which will 
give clear insights into the anatomy of the Higgs-boson properties, the low-energy observables, 
and the DM relic abundance in the region of MSSM parameter space that predicts a signif- 
icantly enhanced diphoton signal. Our analytic analyses will be complemented in Section |3] 
by a detailed numerical study of the light stau scenario, employing state-of-the-art computer 
codes. A summary of our main results and conclusions are presented in Section |4] 



2 Analytic Results 

In this section we present simple formulas that allow to understand the impact of the various 
MSSM parameters on the prediction of the Higgs-boson properties (its mass, production cross 
section, and decay rates) as well as on the expectations for the most important low-energy 
observables {i.e., B — > ^^7, Bg — fi'^fi', and {g — 2)^). The requirements to achieve a proper 
thermal DM relic density are also discussed. Since in our numerical analysis of the MSSM 
parameter space, presented in Section [3| we will employ the most advanced calculations, 
including all relevant contributions, the analytic expressions presented hereafter serve mostly 
an illustrative purpose. 



2.1 Anatomy of Higgs-Boson Mass 



In the decoupling limit of the MSSM, i.e., M\ ^ M|, the lightest CP-even Higgs boson (for 
a review including original references, see [l9]) acquires the squared tree- level mass 



Ml 



Mi c: 



2 „2 



Z 



Mz 2 



(1) 



Here and in what follows we employ the shorthand notations S2/3 = sin (2/3), C2/3 = cos (2/3), 
etc. Because oc 1 in the large-t/? limit and given the minus sign of the decoupling correction 
proportional to M|/M^, it follows that M| < M|. Yet, the scalar field h has SM-like couplings 
when the pseudo-scalar Higgs-boson mass Ma is large, so that this state should have been 
discovered at LEP, if it were not for the radiative corrections which push its mass upward 
from the tree-level upper bound of Mz to Mh > 114.4 GeV 20 . 

In fact, these higher-order corrections can be very large, since the scalar sector of the 
MSSM involves strong couplings, such as those to the top quark and its scalar partners the 
stops. In the limits M^,t^ — )■ 00, that are the relevant ones for the upper bound on Mh, these 
corrections are simple to evaluate. The dominant one-loop contribution to ([T| arises from an 
incomplete cancellation of top-quark and top-squark loops 21 23 , and can be approximated 
by 
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where L^i = In (m^/m?) with m? = m^-^mj^, while = At — fi/t^ denotes the stop-mixing 
parameter, which depends on the trilinear stop-Higgs boson coupling At and the higgsino mass 
parameter fi. We infer that for fixed stop spectrum scale m^-, the Higgs-boson mass correction 
from top/stop loops is maximised for \Xt\ = y/Errit ~ 2.4 m^, which is referred to as "maximal 
mixing" scenario. The manifest symmetry of ^ under the sign flip Xt — )■ —Xt is broken by 



finite two-loop threshold corrections, which induce a term Aas/i^Xt/mi 29,30] that leads to 
slightly larger values of the Higgs-boson mass for XtM^ ^ AtM^ > 0. 

For large tfs there are further contributions to ([T]) that can be relevant [26 -28 . These 
corrections arise from the sbottom and stau sector and take the form (/ = b, f) 

where = 3, A*"^ = 1, = mj_^mj^, and we have ignored logarithmic terms for sim- 
plicity. In the limit of interest, the bottom (tau) Yukawa coupling receives important one- 



loop correction whose dominant contribution depends on sgn (/xMs) (sgn {11M2)) 131 34 . The 



choice /iMs > (//M2 > 0) tends to reduce the tree-level Yukawa coupling yb = y2mh/{vcR) 
{Ur = \/2mr/{vCfj)), and as a result decreases the strictly negative sbottom (stau) effect te| 
on the Higgs-boson mass. 

The relations ([T| to (|3]) allow to draw some general conclusions concerning the impact of 
the observation of a Higgs boson with Mh ~ 125 GeV on the MSSM parameter space. First, 
from the expression M|, we deduce that large values oitp and Ma are crucial to raise the tree- 
level Higgs-boson mass as much as possible. In order to achieve a shift of AM| (85.5 GeV)^, 
the top/stop effects ^ have to be very large, which typically requires > 1 TeV (heavy stop 
spectrum) and/or |^t| ^ 2 TeV (large stop mixing). Since ([3]) scales with the fourth power of 
/i, the value of the higgsino mass parameter should not be too large. Finally, for sgn {AtM^), 
sgn (/iMs), and sgn {11M2) equal to -|-1, subleading negative corrections to M| are minimised. 
We will see below that (some of) the mentioned MSSM parameters also play an important role 
in the production and the decay of the Higgs boson. This fact leads to interesting correlations. 

2.2 Anatomy of Higgs-Boson Production and Decay 

An elegant way of obtaining the interactions of a light Higgs boson, is to construct an effective 
Lagrangian by integrating out heavy degrees of freedom. The resulting effective Higgs-boson 
couplings can be found most easily by utilising low-energy theorems, which relate amplitudes 



with different numbers of zero-momentum Higgs-boson fields pj5p^ (see 37 38 for the original 
idea). Here we will apply this general framework to illustrate the main features of Higgs-boson 
production and decay in the MSSM. In order to present transparent formulas, we will again 



focus on the leading corrections in the decoupling limit 39 . 

In the SM, Higgs-boson production via gluon-gluon fusion receives its dominant contribu- 
tion from triangle diagrams involving top quarks. In the limit of a light Higgs boson, in which 
we are interested in, the corresponding form factor can be replaced by its asymptotic value 

Fi/2(n) ^ lim Fi/2(r,) = 1, (4) 
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where = Am^/M^. In the infinite mass hmit, one has furthermore 19 



Fi/2(r,) 4 ^ ^ 

with -Fo(ri) encoding the effects of scalar loops. 

In the MSSM, the modification of the Higgs-boson production cross section in gg — )■ h, can 
be approximated by 



= 7 7-\ = (1 + S) ~ 1 + / > 



(6) 



where and represents the effects of top-squark and bottom-squark triangles, respectively. 
Notice that ^ ignores the fact that the top-quark and bottom-quark Yukawa couplings in the 
MSSM differ from those in the SM. As will become clear later on, this mismatch is subleading 
in the M|/M^ expansion, and hence can be neglected for our purposes. 

The power of the Higgs low-energy theorems arises from the fact that in the decoupling 
limit the corrections be simply obtained by differentiating the mass-squared matrices 

A1?^ with respect to the mass of the corresponding SM quark. For the top squarks, one has 



where ^rig^^^^ are soft SUSY-breaking masses and Dq^^u^ = O^M^)- A similar expression holds 
in the case of the bottom-squark sector. Ignoring the contributions from the D-terms, which 
are numerically subleading, the master formula for the top-squark contribution to ([6]) reads 

-]^t^ In [det (Ml)] ^ ^ I ^ + ^ ?V I ' (8) 

4 ^ dm; l v tyj 4 i ™ 2 2 ^2 ^2 1 ' v / 

where the multiplicative factor 1/4 in the first line stems from the normalisation ([s]) of the 



scalar form factor, and the final result agrees with the expression given in 40 41 . From (|8j) one 
infers that the amount of mixing in the stop sector, parametrised by Xt, determines whether 
the ratio ^ is smaller or larger than 1. For no mixing, corresponding to Xf = 0, one has 
-R/i > I5 so that Higgs-boson production in gg ^ h is enhanced with respect to the SM. On 
the other hand, if Xt is parametrically larger than the mass eigenvalues m^^ ^ (with m^-^ < m^-^) 
of g, then Rh < 1, meaning that the Higgs boson is less likely to be produced. The fact that 
in the MSSM, in order to make the Higgs boson sufficiently heavy, one needs large/maximal 
mixing, i.e., \Xt\ ~ m^^ m^^ , then tells us that for a random MSSM parameter point that 
gives M/j ^ 125 GeV one should find a suppression of o-{gg h). In fact, this is precisely what 
happens. As a final remark, we add that the sign of the new-physics corrections to Higgs-boson 
production in gluon-gluon fusion is in many models closely related to the (non-) cancellation 
of the quadratic divergence in the Higgs-boson mass |42]. From g and g, we see that in 
the MSSM there is a strong anti-correlation between and Rh, driven by Xt. 
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The calculation of the sbottom contribution to proceeds along the lines of (|8]). Since 
ml <C m^, notable effects can only arise if the mixing in the sbottom sector is very large. In 
this limit, one has approximately 

(9) 



4m? m? 

01 02 



where Xf, = Ab — fitp. Obviously, this correction is strictly destructive and can only be 
important if either the trilinear term At and/or the combination fitp is sufficiently larger than 
the sbottom masses mr . 

"1,2 

In order to describe the decays of the Higgs boson, we define the corrections factors 



T{h VV) 
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T{h ^ VV) 



SM 



r{h ^ //) 
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SM 



(10) 



for V = W,Z, 'J and f = b, r, t. 

With respect to rtsj), the derivation of the leading non-decoupling corrections to Kyj is 
complicated 35 ,36 ,41 by the fact that in the MSSM one has not a single, but two neutral 



scalar fields which develop a vacuum expectation value (VEV) 

1 , . „ 1 



V2 



{Vu + Cah + SaH) 



V2 



Oct 



Here (-7r/2 < a < 0) 



a 



arctan t 



'2/3 



M\ + Ml 



Ml 



(12) 



and h (H) denotes the lighter (heavier) CP-even Higgs-boson mass eigenstate. Furthermore, 



Vu/vd = t/s and ^/vl + 



246 GeV, so that v = Vu/sp = Vd/cp. 



In the case of Ky/, which encodes the modification of the Higgs-boson coupling to a pair of 



W bosons, the presence of the mixing angle a in (11) leads to 



■1 + 
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8M\ 
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(13) 



where we have used that Mw = g/2 \/v\^ v\ and in the last step expanded (12) in powers 
of M\IM\^ retaining only the first non-zero term in the Taylor series. An analogue formula 
applies to k^, which implies that the Higgs-boson couplings to massive gauge-boson pairs 

Kz- This universal correction (13) is strictly 



are affected in a universal way, Ky 
destructive, but of order M\IM\^ and thus numerically insignificant as long as one sticks to 
the decoupling limit. Since s\p oc 1/t^, it is further reduced in the large-t^j limit. 

In the case of the Higgs-boson couplings to tau leptons only the second term of the differ- 
ential operator in (13) contributes because m,- = f^/v^?/^. One obtains 



-1 - 



C/3 



2M| 2 



M\ 



(14) 



The correction is obtained from ( 14 ) by the replacements — Sq, — )■ c^, C/j — s^, and 
Since vfit = Vu/V^Vt only the first term in the bracket of (13) results in a correction to Kt- 
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We see that the shifts in the tree-level couplings of the Higgs-boson to fermion pairs fall off 
quadratically in the limit M\ ^ M|, and that /tr > whereas Kt < 0. In the large-t^ limit 
Kr oc 1 and nt oc 

The coupling of the Higgs boson to bottom quarks receives important t/3-enhanced loop 



corrections involving charginos and gluinos 33,43]. In the limit ^ 1, we find including the 
leading terms 

1 Mi + {AMl), + Ml 

"^^rr^ Ml ' ^^^^ 

where the expression for the tree-level Higgs-boson mass and the dominant one-loop correction 
are given in ([T]) and (|2|, respectively. Furthermore, 



with = m?//i2, = m?/M|, and 

/(^) = - 



fix 



2a,/iM3 




(16) 
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X 
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X 



x] 



(17) 



Notice that f{x) is positive definite with /(O) = 0, /(I) = 1/2, and /(oo) = 1. This feature 



together with the appearance of the combinations iiA^ and i^M^ in (16) makes the correction 



introduced in (15) quite sensitive to the choice of MSSM parameters. While the modifications 
Kb^t have only a minor impact on Higgs-boson production, as already anticipated in (|6]), we 
will see below, that the decoupling corrections to Kb can play an important role for the decays 
of a light Higgs boson. 

In order to give an explicit result for k^, we recall that within the SM, the process h — )■ 
77 is dominated by the virtual loop-exchange of W bosons. These contributions interfere 
destructively with the top-quark amphtude. One has 



Fi{tw) 



lim 



6.24 



(1^ 



where Fiiji) is the form factor associated with vector-boson loops 19 and the numerical result 
corresponds to our reference value = 125 GeV for the Higgs-boson mass. With the help of 
F]Y, we write the modification of the Higgs-boson coupling to two photons as 



1 



1 



- - 



Kr — Kf + Kh± + 



(19) 



where the stop and sbottom contributions, i.e., ki and k^, have already been given in 
and (|9|. Notice that compared to ^ they appear here with opposite signs, signalling that 
constructive interference in gg ^ h goes along with destructive interference in /i — 77 and 
vice versa. This correlation between the loop-induced effective hgg and /177 couplings is a 
general feature in new-physics models with coloured fermionic/scalar partners. 
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The diphoton channel also receives contributions from stau, charged Higgs-boson, and 
chargino loops. The stau corrections takes the form 



Ami rnl 



(20) 



with Xj- = AT-—^tp. Like the correction (20) is only important if the stau-mixing parameter 



satisfies ^ ^ and the lighter stau mass eigenstate is not too heavy. The former 
requirement demands both tp and to be large, and in this region of parameter space, stau 
loops necessarily lead to an enhancement of T{h — )■ 77). 

The charged Higgs-boson effects are, on the other hand, strictly destructive in the MSSM. 
In the decoupling limit, we find 

« [Ml- - -^M14,) m (Ml.) - - . (21) 

where M^± = M\ + M^. Because the spin-zero amplitude Fo(ri) is suppressed by a factor of 
1/4 relative to -Fi/2(Ti) (see ([5|), and M^± k, M\ in the decoupling limit, the correction (21) 
has only a very minor effect on the diphoton decay. Note that cL — )■ 1 for tp ^ 00. 



The last ingredient in ( 19 ) is provided by triangle graphs with internal chargino exchange. 



In terms of the chargino mass matrix 

the corresponding coefficient can be written as 

d 2M^ 
^x^ ^ -M'w In [det (-Mj±-Mx±)] ^ ^gn [det {M^.)] ^ s,, . (23) 

^ Xi X2 

Since for sufficiently large values of M2 and /i, one has det (A^^±) ~ /iM2, the overall sign of 
K^± is determined by the one of the product of the higgsino and wino mass parameters. It 
follows that for /iM2 > (/iM2 < 0), charginos enhance (suppress) the /;,—)■ 77 rate. Since 
S2/3 (X 1/tp the effects are largest for low tp, and numerically very important if the chargino 
spectrum is light. 

When converting the above results into branching ratios, one must bear in mind that the 
total decay rate T{h) of a light SM-like Higgs boson is dominated by its decay into bottom 
quarks. For Mh = 125 GeV, one has BR{h hh) ^ 60%, BR{h WW) ^ 21%, BR{h 



gg) ^ 7%, BR{h r+r") ^ 6%, and BR{h ZZ) ^ 3%. It then follows from (10) that 



i?r = ^ 0.60 (1 + Ktf + 0.07 (1 + Kgf + 0.33 . (24) 

J- WSM 

Note that only the shift Kb has been included here, while the tree-level corrections and kw,z 
have been neglected. This is a very good approximation, since BR{h — )■ bb) ^ BR{h — )■ r+r") 



and (13) is relative to (15) suppressed by an additional power of M^/M\. 
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At this point we are ready to work out the products of the production cross section times 
branching ratios for the various Higgs-boson decay channels. These are the key observables 
that will be affected by the different MSSM contributions. Defining 
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X 



[a{pp h) BR{h X 



MSSM 



1 



i 1 + 1.86 Kg - 1.20 /tb + 2Kx 



(25) 



J] (1 + 

i=g,X 



we find the following semi-analytic results for the most interesting final states X containing 
either massive vector bosons 
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1 + f^htp 
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(26) 



where V = W, Z or diphotons 



R^^l + 0.33 
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0.43 



2 2 
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+ 0.10 
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+ 1.63sgn (/iMa) 



1.20 



Ml + {AMI 



t + Ml 



(27) 



1 + (^btp 



Ml 



The result for the relative signal strength of the hh channel is obtained from (26) by 



simply replacing —1.20 by 0.80. In order to obtain the above expressions we have included 
all non-decoupling corrections, z. e., ([s]), (|9]), (20), and (23), but apart from (15) neglected all 
contributions that vanish in the limit Ml ^ M|. We furthermore took the limit — )■ oo, 
keeping only the leading corrections, and for simplicity replaced sgn [det {Jvl^±)] by sgn (/iM2). 
We also remark that measurements of the double ratio R^/ Rw,z ~ 1 + 2 /t^ would allow for a 
clean extraction of the /i — > 77 amplitude, since it is independent of Hr,b (see 
detailed study). 



44 for a recent 



The formulas (26) and (27) are the main result of this section. They exhibit interesting 
correlations with the expressions presented in Section 2A_ Focusing first on the correction 
{AMDi introduced in ([2|, we see that in the limit of maximal stop mixing \Xt\ ~ \/6m^^7n^ 
(needed to lift the Higgs-boson mass to 125 GeV) the ratios Rw,z and i?^ are necessarily re- 
duced. In terms of m? = m^_^mi^, the shift in Rw,z {R-y 
(— 1.3 m^/m?), which amounts to a correction of around 



is given approximately by —1.9 rrq/rn^ 
1% (-4%) for mi = ITeV. Large 
mixing in the bottom-squark sector will further suppress the latter ratios. The decoupling cor- 
rections affecting the Higgs-boson couplings to the bottom and tau go in the same direction. 
Numerically, one finds a universal shift of —2% for M^ = 1 TeV. Positive corrections to R^ 
can arise from chargino loops if sgn (/iM2) = +1, which, as mentioned below ([3]), helps also to 
diminish the negative correction (AM|)f to the Higgs-boson mass. Since (23) is t/3-suppressed, 
there is however a generic tension between large chargino effects in R^ and saturating the up- 
per limit on the tree-level Higgs-boson mass following from ([T]) . A way to enhance R^ without 
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running into immediate problems at the tree level, is provided by a light stau with large mix- 
ing Xj. K, fitp, which requires that both and fi are large. To give an example, employing 



m| = rrir^mf^ = (200 GeV)^, t/3 = 50, and /i = 1 TeV in (27), one finds that is changed by 
+50%. Since the correction (AM|)^ - in (3) is proportional to — /i'^/m?,, one expects however 
an anti-correlation between the size of trie stau contribution to Rj and the loop-corrected 
Higgs-boson mass. As we will discuss in the next section, some of the MSSM corrections to 
low-energy observables are also changed significantly in the limit — )■ 00, which leads to 
further correlations, strengthening the constraints on the parameter space. 

2.3 Anatomy of Low-Energy Observables 

The discussion in the previous two sections should have made clear that the part of the 
MSSM parameter space with M^^tp — )■ 00 represents a phenomenologically interesting region 
for Higgs-boson physics. In the following, we will extract the dominant MSSM corrections to 
B — )■ Xs7, Bs — yU^/U^, and — in this limiting case, highlighting the existing correlations 
with M/i, Rw,z^ and R^. 

We start our discussion by considering the inclusive radiative B Xg'-y decay. Taking 
into account only the most important corrections due to the Wilson coefficient C7 of the 



electromagnetic dipole operator leads to the following ratio 45 47 



Rxs = Tjp/R ^ Y \ ^ 1 -2.61 AC7 + 1.66 (AC7) , (28) 

where AC7 represents the additive correction appearing in the high-scale Wilson coefficient 
C7 = Cf^ + AC7 ^ -0.19 + AC7 |48). 

Within the MSSM, the Wilson coefficient AC7 receives important contributions from loops 
involving tops and charged Higgs bosons, AC^ , as well as stops and charginos, AC7 . In 
the decoupling limit with ^ 1, the former corrections can be approximated by 



where Ltfj± = In {mf/Mfj±) with M^± ^ M^. Subleading t/j-enhanced corrections dominated 
by gluinos 49| 50 have not been included here. Such effects enhance (suppress) /S.Cj^ for 
negative (positive) values of /iMa, but are irrelevant for our considerations. For = 1 TeV 
one finds ACf ^ ^ -0.03, where the minus sign refiects the well-known fact that in — >■ Xg'y 



the charged Higgs-boson corrections interfere constructively with the SM amplitude 51,52 



Since the charged Higgs-boson effects (29) are not tp enhanced they render only an in- 



significant correction in the MSSM parameter region of interest. The dominant one- loop 



contributions to (28) are therefore provided by diagrams with up-type squarks and charginos. 
These grow linearly with tjj and take the form 

2 

AC^ ^-^A,tp^^g{xi^). (30) 
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Here 



12(1- x) 



6(1 



Inx . 



(31) 



The result (30) agrees with the 
Since the function g{x) is strictly positive for x g]0, oo[ with g{l) = 



Notice that we have included above only the correction due to top squarks and higgsino-like 
charginos, while the wino-like contribution has been suppressed. In the limit of large higgsino 
mass parameters, that we are mainly interested in, this is a good approximation, because 
the latter corrections scale as M2/ fi M^/m^ for 3> M2. 
expression given in |53j 

5/72 0.07, the sign in (30) implies that for fiAt > {fiAt < 0) the branching ratio is larger 
(smaller) than the SM expectation (33}[54] . The dominant t/j-enhanced gluino corrections, not 
shown in (30), again enhance (suppress) the chargino contribution for sgn (/iMs) = — 1 (+1), 
but leave the qualitative dependence of ACj on fiAt unchanged 49 ,50 . For the choices = 
50, m-t = 1.5 TeV, = ITeV, and |A| = 3TeV, one finds AC^ ^ -sgn (/xAt) 0.12, which 



depending on the sign of iJ,At corresponds to an enhancement /suppression of the ratio (28) 
by C(30%). Shifts of this size in BR(_B — )■ ^^7) are detectable given the present theoretical 
calculations and experimental extractions. 

After B — > Xg'y, we now analyse the structure of the MSSM contributions to another 
"standard candle" of quark flavour physics, namely — )■ fi~^fi~. We begin by defining the 
ratio 



R 



cl) 



(32) 



BR(5, i2+i2-)sM 

where Cs and Cp denote the dimensionless Wilson coefficients of the semileptonic scalar and 
pseudo-scalar operators. In the large-t^ regime these coefficients have the most important im- 
pact on R^+^-. The term linear in Cp arises from the interference with the SM contribution 



Cf}'^ ~ —4.2 55 ,56 to the Wilson coefficient of the semileptonic axial- vector operator. Impor- 



tantly, for Cp > it interferes destructively with the term proportional to (C| + Cp). This 
implies, on the one hand, that for positive values of Cp the stringent bounds on Bs — > 
are more easily evaded and, on the other hand, that a pseudo-scalar contribution of the correct 
size will lead to a suppression of the purely leptonic decay mode below its SM value. 

Within the MSSM the contributions to Cs and Cp with the strongest dependence arise 
from neutral Higgs double penguins 57 . In the decoupling limit, one has Cs ~ — Cp with 



Cp ^ ixAt 



2 M\ 



(33) 



Henceforth we use the shorthand notation sw = sin6w etc. to indicate trigonometric func- 
tions of the weak mixing angle. Notice that our definition of the semileptonic scalar and 
pseudo-scalar Wilson coefficients differs from that of |58| by a factor of rrib. The parameter Sb 



introduced in (16), parametrises loop-induced non-holomorphic terms that receive their dom- 
inate contributions from higgsino and gluino exchange. The Wilson coefficients C5 and Cp 

but these are of no concern as 



58,59 



also receive various other contributions in the MSSM 
long as one is interested in the general structure of the effects only. Recall that /(x) > 0, so 
that the sign of Cs [Cp) is opposite to (follows) that of the combination fiAt. Notice finally 
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that ( [l6| ) introduced a dependence on sgn (yuMa) 

/iMg > (/iMg < 0). 



and that (33) is suppressed (enhanced) for 



From what has been said below (32 ), it should be clear, that the latter feature has important 
implications. In particular, it follows that finding BR(i?s — )■ yU+yU^) close to its SM value leads 
to a two-sided bound on the product fiAt, if double Higgs-penguin contributions provide the 



dominant new-physics effect in the purely leptonic decay. For example, consider R 



Combining (32) and (33), translates into 



0.16 

TeV^ 



< 



{l + e,tpfm\Ml 



< 



1.37 
TeV^ 



< 1.3. 



(34) 



For simplicity we have set x^-^ = 1 here, which corresponds to the case of degenerate masses. 
Clearly, this double inequality provides a non-trivial constraint on the MSSM parameter 
space with ^ 1. Realize that the apparent asymmetry in (34) arises from the fact 
that Cp > is preferred over Cp < 0, because a small positive pseudo-scalar contribu- 
tion leads to a cancellation between the linear and quadratic terms in (32), and that values 
1/(1 + ebtpffiAt/{mlMl) {tfs/bOf ^ 0.6/TeV^ correspond to a reduction of BR(5, /i+/i") 
relative to the SM by about 50%. 

We now turn our attention to the anomalous magnetic moment = [g — 2)^/2 of the 
muon. In the MSSM, this observable receives one-loop contributions from Feynman diagrams 
with neutralino and smuon as well as chargino and sneutrino exchange (see 60,61 for topical 



reviews). The latter diagrams dominate in almost the entire parameter space [62] and can be 
approximated by 



with X 



2 



/ml and 



m 



X 



X 



X21' 



h{x) 



X — 3 
(1^ 



'1 



Inx - 



X) 



(35) 



(36) 



^2jy by convention and h{x) is a monotonically decreasing function of x with 
0, it follows that the sign of (35) is determined by that of fiM2. 
: |y(i| = ITeV, IM2I 



< X 

h{l) = 2/3 and h{oo) 
Employing = 50, rrii, 

finds a shift Aa};^ ^ sgn (yuMg) 2.2 



300 GeV, m + ~ IM2I, and m ± ^ lul, one 

Xi A.2 

which for yU.M2 > removes almost entirely the 



well-known tension between the experimental result and the SM prediction for a. 



The formulas (30), (33), and (35) provide a good starting point to discuss the correlations 



between the low-energy and Higgs-boson observables, which we expect to see below in our 
numerical analysis. Choosing, for definiteness, sgn (Mj) = +1 with i = 1,2, 3, we infer from 
Sections 2.1 and 2.2 that the parameter region with large t/s and large and positive At and fi is 



well suited, on the one hand, to explain the relative heaviness of the Higgs boson, and, on the 



other hand, to allow for a significant enhancement of the /i — )■ 77 signal. From (30) it follows 



that the above parameter choices inevitably lead to an enhancement of the branching ratio 
of 5 — >^ -^^57. While the Bg — )■ fi'^n~ decay rate can in principle be both enhanced as well as 
reduced for /lAt > 0, we will see in the next section, that for yuMa > one typically observes 
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a suppression of the Bg — )■ fi^fi branching ratio whenever the diphoton signal is enhanced. 
Notice that the presence of the t/j-enhanced corrections in (33) plays a crucial role in this 



context. In fact, from (|34|) it is evident that parameter choices that give e;, > are favoured 

+1 in the considered parameter region, 
data, if the sneutrinos are 



over those that lead to < 0. Finally, since sgn (jj,) 
one also expects a visible improvement in the description of the 
sufficiently light. 



2.4 Anatomy of DM Relic Abundance 

Another appealing feature of the special MSSM parameters under consideration, is the possi- 
bility to generate the correct DM relic density 

^ 1.07 ■ 10^ xj 

^ GeV MpiV^a., ' ^'^^ 



through stau coannihilation 63 66 with the lightest SUSY particle (LSP), represented by the 
lightest neutralino mass eigenstate. In the above equation, Mp\ ^ 1.22- 10^^ GeV is the Planck 
mass, ~ 81 is the number of relativistic degrees of freedom at the freeze-out temperature 
Tf, and Xf = m^o/Tf ^ 22 for the masses of the lightest neutralino we are interested in. 
In the case at hand, the effective cross section 

(Tofr ^ "xx'^xx + "x^«x^ + ^ (1/3^/) ' (38) 



entering (37) receives two important contributions. The first one arises from neutralino anni- 
hilation into tau pairs, induced by stau exchange in the t and u channel. For a mostly bino-like 
neutralino, m^o ^ |Mi|, and large mixing in the stau sector, we find the expression 

1 

''^^'^32^^^'^''"^^' ^^^^ 



where r^^ = mf^/m^o and 

^■('') = JT^. ■ («) 

In order to achieve a pronounced h ^ '-yj signal, the stau has to be light and strongly 
coupled to the Higgs. It turns out that in this case the dominant coannihilation channel is 
XiTi — )■ /ir, which proceeds through the t-channel exchange of the lighter stau state jij. Its 



contribution to (|38|) is approximately given by 



a 



X 



with r^x = mh/^Xi 



k(run)^ (r. - r, + + „ + IP^ ^ ^^^^ 

ri(ri + 1) (r| - ri(ri + 1)^) 
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The results (39) and (41) can be derived from the general expressions given in \66\. Notice 
that the coannihilation contribution a^f is proportional to m'^X^. Like in (20) this factor 



arises from the fifih coupling. 

The last missing ingredient needed to allow for a qualitative understanding of flBuh"^, 



the Boltzmann weight factors a^^ and a^f appearing in (38). In the parameter space that 
allows to reproduce the observed relic density, we obtain 



a^^^l, a^,^-e2°-^(i-^x). (43) 



The expressions (43 ) show that the size of coannihilation is exponentially sensitive to the mass 
splitting rrir^ — fri^o, meaning that the process Xi^i hr only becomes important in regions 
of parameter space that permit values of Vf^ close to 1. 



Based on (37) to (43), we are now able to discuss the restrictions that the observed relic 
density fl-Duh"^ ~ 0.11 imposes on the MSSM parameter space. For neutralino masses m^o < 
100 GeV, it turns out that XiXi ^ t^t^ represents the dominant DM annihilation channel. 
In this case, one has f^DM^^ ~ (^dm^^)xx with 

{n^Mh\,^iA.io~' {^^)\^ + 4,r- (44) 

It follows that in order to achieve a large enough DM abundance one needs a sizeable mass 
splitting rrif^ — m^o. For example, for m^o ^ 30 GeV masses of ruf^ ^ 90 GeV result in viable 
^DM^^ values. Increasing the neutralino mass requires slightly larger splittings m^^ — m^o 
that reach a maximum at m^o ~ 60 GeV. At this point the mass difference starts decreasing, 
leading to ruf^ — m^o fa 35 GeV for m^o ^ 100 GeV. 

For somewhat heavier neutralinos m^o > 100 GeV the proper DM density requires mass 
differences m^^ — m^o of order 20 GeV, and as a result (jllj) has to be included. In consequence, 
^DM^^ ~ (f^DM^^)xx + (^DM^^)fx "with the latter contribution approximated by 



inr,uh%, ^ -2.5 e--(--^) . (45) 

Notice that the exponential suppression of (fiDM^^)fx ^"7 f^fi — '^^ specific MSSM 

scenarios we consider balanced by the large mixing parameter Xj. ^tp. 



3 Numerical Analysis 



We now turn to the numerical analysis of the Higgs-boson and the low-energy observables as 
well as the DM relic density. All results presented below have been obtained using Soft- 
SUSY 3.3.3 [67] for the spectrum calculation, HiGLU 3.11 l68] and Hdecay 4.45 feF] for the 



computation of Higgs-boson production and decay, and SuperIso Relic 3.3 70|71j for the 
calculation of the relevant low-energy observables and the DM abundance. We have compared 
our results against SuSpect 2.41 [72] and FeynHiggs 2.8.6 [73], and found good overall 
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Figure 1: Results for Mh (upper left), Rb (upper center), (upper right), R^ (center 
left), nif^ (center), Rx^ (center right), R^+^- (lower left), Aa^ (lower center), and 
^^DM^^ (lower right) in the = 60 scenario. The Higgs-boson and lighter stau masses 
are given in units of GeV. The dotted black lines indicate the parameter regions with 
R^ > 1, while the dashed black lines correspond to the 95% CL regions favoured by 
B — )■ Xg'j. See text for details. 



agreement between the different programs for most of the observables. The biggest numerical 
differences arise for the Higgs-boson mass and its branching ratios. In the former (latter) case 
we find relative changes of typically below 5% (10%). The observed differences can be traced 
back, on the one hand, to the use of different renormalisation prescriptions (DR vs. on-shell 
scheme) and, on the other hand, to the different treatment of higher-order perturbative cor- 
rections. The quoted relative errors give an indication of the theoretical uncertainty plaguing 
our calculations of the Higgs-boson observables, and we will comment on its impact on our 
numerical analysis below. For a detailed comparisons between the publicly available programs 
dealing with the mass of the Higgs boson in the MSSM, we refer to |74J . 
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To begin with, we focus on scenarios with = 60. The other choices of the MSSM 
parameters are Ma = ITeV, Mi = 50GeV, M2 = 300 GeV, M3 = 1.2 TeV, rhg^ = mu-, = 
1.5 TeV, rfiL^ = fhi^ = 350 GeV, while we take common soft SUSY-breaking masses of 1.5 TeV 
and 2 TeV (ITeV) for the remaining "left-handed" and "right-handed" squark (sleptons). We 
furthermore employ Ab = 2.5 TeV and Ar = 500 GeV, while the first and second generation 
trilinear couplings take the same values as those of the third generation, as they have essentially 
no impact on the observables in question. Hereafter we will refer to this specific choice of 
parameters as the "t/j = 60 scenario". In Figure [l] we show the results of our numerical scans 
in the Af-fi plane. We restrict ourselves to the quadrant with At > and /i > 0, which 
shows the most interesting effects and correlations and is the only one that allows for a good 
description of all data. 

From the prediction for the Higgs-boson mass, we see that for the above choice of SUSY 
parameters, the trilinear term At has to lie in the range of [2, 5] TeV in order to push up to 
[123, 129] GeV. The latter range is allowed by the ATLAS and CMS data (l]|2], if one accounts 
for the parametric errors from the SM input (with the largest uncertainty arising from the 



top-quark mass, mf"^*^ = (173.3 ± 2.8) GeV 75 ) as well as the theoretical uncertainties in the 



MSSM calculation of the mass of the Higgs boson 74 . Of course, the need for large trilinear 



stop-Higgs boson couplings is an immediate consequence of ([T| and p]). The anti-correlation 
between the Higgs-boson mass and the fj. parameter, as implied by ([3|, is also clearly visible 
in the panel. We emphasise that the shown predictions correspond to DR input parameters 
and that the results obtained in the on-shell scheme (as used for example in FeynHiggs) 
would differ to some extent. In particular, the values for At needed to accommodate a Higgs- 
boson mass consistent with the LHC observations can be smaller by a factor of up to 2. This 
pronounced dependence on the renormalisation prescription should be kept in mind when 
interpreting our numerical results. 

We now turn our attention to Higgs-boson production and decays. Since we are in (or close 
to) the maximal stop-mixing regime, we expect from ([s]) that the Higgs-boson production cross 
section should be suppressed with respect to the SM. In fact, we obtain ~ 0.95 (and stop 
masses m^^ ~ 1.3 TeV and rrit^ ~ 1.6 TeV) throughout the entire parameter space depicted 
in the panels. The deviations in the relative signal strength Rb of the Higgs-boson decay to 
bottom quarks are more pronounced than those in Rh and correspond to enhancements of 
roughly 20%. Notice that these shifts are due to the terms and (AM^)^ in (15), the latter 
of which introduces a positive correlation between Rb and M^. This feature clearly manifests 
itself in the two panels. As a consequence of the suppression of Rh and the enhancement 
of Rb, the Higgs-boson decays to W-, Z-boson, and tau pairs are all suppressed. For the 
considered MSSM parameters, we find Rw,z ~ 0.7 and Rr ~ [0.65,0.80]. We add that Rr 
shows a noticeable dependence on the higgsino mass parameter that stems from chargino and 
neutralino effects in the tau Yukawa coupling. While the structure of these corrections is 
similar to those in (16), these contributions have, for simplicity, not been included in the 
approximate formula (14). 



As anticipated after (20) and (27), the prediction for the diphoton signal R^ depends very 
strongly on the amount of mixing Xr ~ fitp in the stau sector. We observe that in the 
studied scenario, values of i?^ > 1 can only be obtained in a narrow mass window around 
II ~ 900 GeV. For our choice of soft SUSY-breaking masses tHl^ = rhi^ = 350 GeV, such large 
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Figure 2: Left: Average /;,—)■ 77 partial decay width normalised to its SM value as a 
function of the lighter stau mass. The shown vertical bars indicate the one standard 
deviations from the averaged values. Right: Values of the lighter stau mass necessary 
to obtain the correct DM abundance as a function of the lightest neutralino mass. In 
the parameter space below (above) the blue region the predicted values of f^oM^^ are 
below (above) the 3.5cr WMAP bound. 



H parameters lead to values of m^^ that are close to the LEP bound of rrif^ > 81.9 GeV 76 . To 



further illustrate the latter feature we display in the left panel of Figure [2] the maximal value 
for the /i — )■ 77 decay width normalised to its SM value, that can be obtained for fixed m^^. 
In order to arrive at the plot we have used the results from a general 19 parameter scan in the 
phenomeno logical MSSM, as described in 77,78 . The fast decoupling of the stau corrections 
in r(/i — )■ 77) is evident from the figure. 

We have seen that achieving a pronounced /i — ?■ 77 signal requires the presence of a 
large fifih coupling. Such an interaction can be potentially problematic, since it may trigger 
additional minima in the scalar potential, and as a result the electroweak-breaking vacuum 
can become metastable 79 82 . A tree-level analysis of the stability of the vacuum, taking 
into account large left-right mixing in the stau sector, leads to the following constraint ISO] 



\fi\t0 < 38.5 (v^+ V^)' - lOTeV. 



(46) 



Although this result cannot be fully trusted as loop effects are very important in the specific 
MSSM scenarios considered in our work, it is clear from Figure [T] that in the slice of parameter 
space that features significant enhancements of the diphoton signal, the bound (46) is violated 
by around 30%. In order to understand to which extend vacuum stability considerations 
exclude a light stau explanation of i?^ ~ 1.7, a full one-loop analysis of stau effects in the 
scalar potential would be required. Such a study is beyond the scope of this article. 

We start the discussion of our numerical results for the low-energy observables with B — )■ 
X,7. Adding to the uncertainty of the SM prediction BR{B X,7)sm = (3.08 ± 0.24) • 
10~^ [45 ,46 ,70] an intrinsic MSSM error of 0.10 as well as the error of the experimental world 

(3.43 ± 0.22) 



average BR{B Xs7)cxp 
level (CL) bound 



10"* 83 



Rx^ = 1.11 ±0.11, 



leads to the following 68% confidence 

(47) 
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where the individual uncertainties have been added in quadrature. At 95% CL one obtains 
instead Rx^ = [0.89,1.33]. The border of this 95% probabihty region is indicated in the 
panels of Figure [l] by the dashed black curves and arrows. From the panel showing Rx^ , 
one infers that for the depicted parameter choices, BR(i? — t- Xs7) is always enhanced (by 



about 20% to 60%) with respect to the SM expectation. As expected from (30), the MSSM 
corrections grow with At and become too large for > 4.5 TeV to allow for an agreement with 
BR(i? — )■ Xs'-j) at the 95% CL. This is an interesting and potentially important finding, since 
the parameter space disfavoured by i? — partially overlaps with that preferred by other 
observations/measurements. In fact, we see from the figure that the B — > constraint 
starts cutting into the already narrow regions in the At-^ plane with Mh ~ 125 GeV and 
R^ > 1. 

A glimpse at the predictions for R^+^- also shows that there is an intriguing correlation 
between BR(i?s — fj^~^fi~) and BR(i? — Xsj). Notice first that the obtained branching 
fractions of Bg — ^ fi~^fi~ are all fully compatible with the bound 

R^+^~<1.5, (48) 



which derives from the 95% CL exclusion BR{Bs — )• /i^/i~)exp < 4.2-10"^ 84 and the untagged 
SM branching fraction BR{Bs — )■ /U"'"/i~)sM = 3.9 ■ 10~^ [70,85,86] (corresponding to a CP- 
averaged branching ratio of 3.5 ■ 10~^), after including a theoretical 68% CL error of 20% 
that is thought to cover both the SM and SUSY uncertainties. Second, we observe that 
essentially all solutions that satisfy B — )■ Xs'~f feature a suppression of Bs — > yu^/^"- In fact, 
asking for both an agreement with Rx^ at the 95% CL as well as an enhanced diphoton rate. 



implies -R^+^- ~ [0.6,1.0]. Notice that given the 1/M^ dependence of (32) and (33) the ratio 
Rfi+^i- is a very sensitive measure of the masses of the heavy Higgses. It follows that the 
deviations quoted above can be reduced by choosing Ma 3> 1 TeV. An observation of the 
purely leptonic Bg decay at the SM level (which seems possible with 2012 LHC data), may 
hence give important insights both on the nature of the /i — )■ 77 excess in the context of the 
MSSM as well as the size of the decoupling scale Ma- 

The last remaining low-energy observable in our study is the anomalous magnetic moment 
of the muon. Despite many changes and improvements in the recent history, the discrepancy 
seen in seems to persist. Combining the experimental result [87j with the SM calculation, 
based on an update of the hadronic vacuum polarisation contributions [88] and the complete 
tenth-order QED corrections [89], results in 

Aa^ = al^P - af^ = (2.5 ± 0.9) ■ 10"^ . (49) 

We see that for the sneutrino masses m^; ^ 1 TeV present in our scenario, the predicted shifts 
in amount to about [1.0, 1.6] ■ 10~^, which leads to a significant reduction of the above 



tension. The observed anti-correlation between Aa„ and fx is readily understood from (35). 



For /i, m[, ^ M2, the chargino-sneutrino corrections to the muon anomalous magnetic moment 



scale as Aa^ oc M2/ fj, M^/mj,. This relation also implies that a notable improvement in (49) 
requires a relatively light slepton spectrum below a TeV, which in our scenario is present due to 
a suitable choice of parameters. The observed correlation between R^ > 1 and Aa^ = O(10~^) 
should therefore not be regarded as a sohd prediction in the entire MSSM parameter space. 
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Figure 3: Predictions for Mh (upper left), Rb (upper center), (upper right), Rx^ 
(lower left), R^+^- (lower center), and f^DM^^ (lower right) in the tp = 30 scenario. 
The mass of the Higgs boson is given in units of GeV. The dotted black lines indicate 
the parameter regions with i?^ > 1, while the dashed black lines correspond to the 
95% CL regions favoured by -B — )■ X^'y. See text for further explanations. 



We add that if slepton mass universality is assumed |10] , the correlation between R^ and Aa^ 
becomes however quite robust. 

Let us now switch gear again and finally examine the predictions for the DM relic abun- 
dance. We see that the obtained values for f^DM^^ range over three orders of magnitude, but 
that agreement with the tight WMAP 3.5cr bound l90] 



f]j3M/i^ = [0.068,0.155] , (50) 



that includes theoretical uncertainties (see for example [91] and references therein), can be 
achieved. In fact, requiring only that the LSP does not overpopulate the universe, i.e., 
flBuh"^ < 0.155, singles out a parameter region in the Af-fi plane that overlaps with that featur- 
ing R^ > 1. The strong anti-correlation (positive correlation) between the ^dm^^ and /i {rrifj 
is also clearly visible in the panels. It is easy to understand these two features by considering 



the pure annihilation contribution (44) to I^dm^^ that effectively limits the size of the mass 
splitting —m^o. Numerically, we find that for our choice Mi = 50 GeV ~ ^x?' require- 
ment of an electrically neutral LSP with ^dm^^ < 0.155 is only fulfilled if ruf^ ~ [80, 120] GeV. 
Since = fhi^ = 350 GeV, such relatively light staus can however only be obtained for large 
/i parameters. The strong correlation between the mass of the LSP and the lighter stau is 
illustrated in the right panel of Figure [2| which displays the parameter region in the m^o-rrif^ 
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Figure 4: Predictions for Mh (upper left), Rb (upper center), Rr (upper right), R^ 
(lower left)), Rx^ (lower center), and R^+^- (lower right) in the = 60 scenario for 
At = 3 TeV. The mass of the Higgs boson is given in units of GeV. The dotted black 
lines indicate the parameter regions with R^ > 1, while the dashed black lines corre- 
spond to the 95% CL regions favoured by S — )■ Xg'j. See text for further explanations. 



plane that is consistent with (50). The shown predictions correspond to the = 60 scenario 
parameter choices with the value of Mi varied. From the figure it is clear that for a fixed value 
of only a narrow range of m^o values is consistent with the WMAP bound. This is turn 
implies that in the light stau scenario a confirmation of the excess in — )■ /z. — )■ 77 will have 
implication for direct and indirect DM searches, since the LSP mass is not a free parameter, 
but fixed to some degree. 

The above discussion should have made clear that a confirmation of the results on the 
Higgs-boson couplings may point towards rather peculiar (and technically "unnatural") MSSM 
parameters, namely large (and positive) values of At, and /i. These special parameter 
choices lead in turn to interesting and testable correlations between various observables. In 
the following, we would like to address the question of how robust these correlations are against 
the variations of some of other MSSM parameters that have been kept fixed so far. 

We start our discussion by studying the impact of t/s. In what follows we employ = 30, 
7711,3 — 170 GeV, and m/3 = 350 GeV ("t^ = 30 scenario"). Notice that the change of the 
soft SUSY-breaking masses is required in order to obtain a very light fi eigenstate, which in 
turn results in a notable shift in /i — )■ 77. Furthermore, this choice of parameters satisfies the 
vacuum stability bound (46). The results of our numerical scans are collected in Figure [3j Let 
us first consider the Higgs-boson mass as well as the h ^ bb and /i — )■ 77 decay signals. We see 
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Figure 5: Predictions for (left), Rx^ (center), and R^+^- (right) in the = 60 
scenario, employing M3 = —1.2 TeV. The mass of the Higgs boson is given in units of 
GeV. The dotted black lines indicate the parameter regions with > 1. See text for 
additional details. 



that while the general pattern of the predictions resembles that obtained for t/j = 60, certain 
differences are clearly visible. First, now even values of /i > 1 TeV lead to allowed Higgs-boson 
masses in the range [123, 129] GeV. This feature is related to the negative corrections (|3| to the 
Higgs-boson mass that scale as (AM|)j^- oc t^. Second, the enhancements in Rf, are slightly 
smaller than those found for = 60, but still amount to shifts in the ballpark of 20%. In 
turn, Rw,z and Rr turn out to be somewhat larger than in the t/j = 60 scenario. Third, the 
enhancements in R^ are significantly smaller now and limited to 15%. This is expected since 
the stau corrections to the effective /i — )■ 77 vertex (20) scale like oc t^. 

Turning to the predictions for B — )• Xgj, Bg -^fj^fi~, and I^dm^^, we first observe that 
in the tfs = 30 scenario, SUSY effects still tend to enhance the rate of the inclusive radiative 
B decay, but that the corrections are about a factor of 2 smaller than before, i.e., Rx^ ~ 
[1.1, 1.4]. Obviously, this is a result of the linear tp scaling of the chargino corrections (30). 
From (35) it follows that the same kind of depletion is also present in Aa^. In the case of 
Bg l^'^ iJ^~ 1 we find that in the parameter space favoured by i? — Xs7, the branching ratio 
of the purely leptonic Bs decay is always suppressed with respect to the SM expectation. One 
has Rfj.+^- ~ [0.6, 0.8]. This feature can again be understood from the interplay of t/j-enhanced 
terms in (33). As in the case of = 60, we finally see that acceptable values of the DM relic 
density can be achieved in the parts of the plane that also give R^ > 1. 

Another important SUSY parameter is the gluino mass M3. In order to study its impact, we 
perform scans in the fi-M^ plane, fixing the value of the trilinear stop-Higgs boson coupling to 
At = 3 TeV and employing the parameters of the = 60 scenario discussed before. Our most 
important findings are illustrated in the panels of Figure |4j In the case of the Higgs-boson mass 
we observe that for fixed fi parameter the predictions for reach a maximum for M3 ^ 1 TeV 
and then start decreasing for increasing gluino masses. This effect is associated to the one-loop 
gluino corrections (16) to the bottom Yukawa coupling that lead to a negative shift in Mh 
via ([3]). In fact, we infer from the predictions for Rj that the most pronounced enhancement 
in this observable occur in a thin stripe with /i ^ 900 GeV and M3 ^ [0.5, 1.5] TeV. Requiring 
Mh > 123 GeV and R^ > 1 hence effectively sets an upper limit on the gluino mass. In the slice 
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of the /i-Ms plane that leads to an enhanced /i — )■ 77 signal, we also see that the predictions 
for h ^ bb, h ^ t^t~, Rxsi ^^id R^+^- vary only moderately with M3. Numerically, we 
find that ^ [1.15,1.21], R^ ^ [0.66,0.76], Rx, ~ [1-2,1.3], and ^ [0.6,1.0]. Notice 

that Rb increases with M3, while Rr and R^+fj,- both decrease with increasing gluino mass. 
These behaviours can be traced back to the t^-enhanced gluino corrections entering (16). The 
observables not explicitly displayed in the figure show either essentially no {Rh, Aa^, and 
^Duh"^) or only a minor dependence {Rw,z) on the gluino mass. 

So far we have restricted ourselves to the case of a positive gluino mass term M3. In 
order to see whether this sign choice has some impact on the obtained results, we repeat our 
numerical scans in the = 60 scenario, using M3 = — 1.2TeV instead of M3 = 1.2 TeV. The 
corresponding plots can be found in Figure |5} From the predictions for M^, we see that for 
M3 < 0, Higgs-boson masses above 127 GeV cannot be achieved. This feature is related to 
the fact that for fiM^ < the tree-level bottom and tau Yukawa couplings are enhanced with 
respect to the case of fxMs > 0, which leads to larger negative corrections (|3| to Mh. We 
add that the Higgs-boson decays to massive gauge bosons and fermions all remain essentially 
SM-like in the studied scenario, while R^ can be enhanced by up to 35%. For what concerns 
the S-physics observables, we observe that fiipping the sign of M3 while leaving the remaining 
MSSM parameters untouched, enhances both BR(i? — )■ Xg'y) and BR(i?s — )■ fi~^fi~). Notice 
that these enhancements originate from the t/3-enhanced gluino corrections appearing in (16). 
While in the case of i? — )■ ^^7 the predicted values of the branching ratio are larger by around 
a factor of 1.1, we find that the branching fraction of Bs — )■ fi~^fi~ is increased by factors 2.5 
to 5. Given the stringent bounds (47) and (48), scenarios with M3 < and large tp, At > 0, 
and n > are hence disfavoured. The remaining observables do not significantly depend on 
the sign chosen for M3, and we hence do not show the corresponding predictions in the figure. 

From the above explorations in the MSSM parameter space, we conclude that many of the 
found correlations are robust, as long as one restricts oneself to the region with > 50, large 
and positive At and fi, ^ 125 GeV, and an enhanced — )■ 77 rate. Positive values of M3 
are also clearly favoured over negative gluino mass parameters. 

Let us finally add that a strong enhancement of R^ can also be achieved by suppressing 
the partial decay rate of the Higgs boson to bottom pairs. While this can be easily achieved 
in the MSSM, a suppression of Rb leads typically to enhanced h — WW, ZZ rates. Given 
that the ratios Rw,z appear to be SM-hke [T||2] an explanation of R^ ~ 1.7 via a suppressed 
h ^ bb width is (at present) not favoured by experiment. 



4 Conclusions 

The announcements of the discovery of a bosonic state by the LHC high-p-r experiments mark 
the beginning of a new chapter in particle physics. While the significance of the various 
measurements is not yet sufficient to tell if the properties of the observed particle agree with 
that of the SM Higgs scalar, the preliminary findings of an enhanced /;,—)■ 77 rate have 
triggered a lot of excitement, in particular, in the theoretical community. With ATLAS 
and CMS accumulating more data, the question of whether new physics or just a statistical 
fiuctuation is responsible for the observed deviation, may be answered by the end of this year. 
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In this article we have studied in detail under which circumstances MSSM scenarios with 
a light stau can give rise to a significant enhancement of the diphoton signal without violating 
other existing constraints from B physics, {g — 2)^, and dark matter. We found that the 
observation of a Higgs-like state with a mass of around 125 GeV combined with the preliminary 
measurements of the Higgs-boson couplings points towards a distinct (but unnatural) choice 
of parameters, namely large values of tan/3. Ma, At, and /x with sgn (A^) = sgn (/i) = +1. 
In this region of parameter space the correct thermal relic density can be achieved, but only 
if one assumes the hierarchy \Mi\ <^ IM2I <^ A typical MSSM spectrum leading to a 
significantly enhanced h '-ff rate as well as the correct value of ^dm^^, therefore contains 
a light bino as the dark matter candidate, a light and maximally mixed stau (often causing 
the vacuum to become metastable), and a heavy higgsino. We also showed that spectra where 
the gluino is much heavier than the squarks can be problematic, since in such a case and 
in turn the enhancements in /i — )■ 77 are bounded from above. 

In the corner of phase space singled out by the Higgs-boson mass, the diphoton rate, and 
the relic density, we found that the predictions for the remaining Higgs-boson observables 
are relatively robust against variations of the other MSSM parameters. In fact, intriguing 
patterns of deviations surface. While Higgs-boson production is typically shghtly suppressed 
with respect to the SM, the h ^ bb rate is generically enhanced by around 20%, which in turn 
results in suppressions of the decays to H^-boson, Z-boson, and tau pairs by a comparable 
amount. Such shifts are in good agreement with the tentative findings by the LHC and 
Tevatron collaborations on the search for the SM Higgs scalar. 

A further consequence of the large and positive values of tan /3, M^, At, and fi are B — )■ Xgj 
branching fractions that are above the SM expectation by about 30%. In view of the ongoing 
effort to improve the theoretical understanding of the inclusive radiative B decay, deviations 
of such an amount may provide a smoking gun signal of the light stau scenario in the future. 
Let us add that in contrast to -B — )■ Xg'y, the predictions for B — )■ Xsl^l~ , B — )■ K*l^l~ , and 
B ^ TV all turn out to be SM-like and well in agreement with experimental results. Similarly, 
the values of the B — )■ D^*^ti' branching ratios are very close to the SM expectations. The 



anomalies seen in the latter channels 92 can hence not be accommodated in the light stau 
scenarios (nor in the full MSSM). We also observe that in the region of parameter space 
favoured by i? — )■ Xg'j the rate of Bg — )■ fJ''^fi~ tends to be smaller than the SM prediction. 
The sign of the gluino mass parameter plays a crucial role in obtaining viable predictions for 
the observables in the 5-meson sector and we found that positive values of M3 are clearly 
favoured over parameter choices with M3 < 0. Since the corrections in Bg — t- ^~ can reach 
up to 40% (for not too large values of Ma), precision measurements of the purely leptonic 
Bg decay, now under way at ATLAS, CMS, and LHCb, might shed further light on whether 
the observed enhancement in — >■ 77 is due to a light stau. We finally showed that under 
the assumption of light soft-SUSY breaking slepton masses, the long-standing discrepancy in 
{g — 2)^ is significantly reduced in the parameter region selected by the Higgs-boson mass and 
the enhanced diphoton signal. 

The choice of the scenarios studied in this paper allowed us to identify the correlations 
which would have been otherwise difficult to highlight in a general MSSM scan. These correla- 
tions might become very useful in case the ATLAS and CMS experiments continue to measure 
a notable enhancement of /i — 77 and/or start to see the first supersymmetric partners. 
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